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We show that electrodynamic dipolar interactions, responsible for long-range fluctuations in mat- 
ter, play a significant role in the stability of molecular crystals. Density functional theory calcu- 
lations with van der Waals interactions determined from a semilocal "atom-in-a-molecule" model 
result in a large overestimation of the dielectric constants and sublimation enthalpies for polyacene 
crystals from naphthalene to pentacene, whereas an accurate treatment of non-local electrodynamic 
response leads to an agreement with the measured values for both quantities. Our findings suggest 
that collective response effects play a substantial role not only for optical excitations, but also for 
cohesive properties of non-covalently bound molecular crystals. 



Polyacene molecular crystals form a fundamental class 
of aromatic solids, and have been extensively studied 
as potential materials for organic electronics [TH3]- It 
is understood that the optical properties of polyacenes 
are very sensitive to long-range intra- and intermolcc- 
ular electrodynamic interactions. This is reflected by 
shifts in the optical absorption frequencies upon increas- 
ing the molecule size or upon solid formation [3], and 
is further exhibited by the visible color of oligoacene 
crystals, which changes from transparent in naphthalene 
and anthracene, to bright orange in tetracene, and deep 
blue in pentacene [4j [5]. The optical absorption spec- 
trum is directly related to the polarizability through the 
Kramcrs-Kronig transformation [B]. Therefore, the ob- 
served changes in the optical spectrum upon crystalliza- 
tion of polyacenes are accompanied by a change in the 
molecular polarizability. In addition, these changes in 
polarization should directly impact the crystal lattice en- 
ergy. However, the effect of electrodynamic intermolec- 
ular interactions on the cohesive properties of molecular 
crystals remains poorly understood. In this Letter, we 
show that the dipolar electrodynamic coupling between 
polyacene molecules reduces the solid dielectric constant 
by 15%, and has an impact of up to 0.5 eV per molecule 
on the computed van der Waals energies and sublima- 
tion enthalpies of these molecular crystals. Our results 
imply that electrodynamic response is crucial for describ- 
ing both the cohesive energy and the optical properties 
of molecular crystals. 

Polyacene crystals are extended aromatic networks 
characterized by polarizable n clouds. Therefore, an ap- 
preciable part of the crystal lattice energy stems from 
ubiquitous attractive vdW dispersion interactions. When 
studying the cohesion of molecular systems, for example 
using density-functional theory (DFT) [7J H] or classical 
potentials [9], the vdW energy is typically computed us- 
ing effective polarizabilities for hybridized "atoms" inside 
a molecule. It is common to approximate the frequency- 



dependent polarizability of every atom using a single ef- 
fective excitation frequency (also called the Unsold ap- 
proximation [10] ) . In this model, the dipole polarizability 
for atom p is written as 



a p [n(r)} 
pV ' l + ( U / Up [n(rW 



(1) 



where a p [n(r)] is the static polarizability of an atom 
p, and u; p [n(r)] is the corresponding characteristic ex- 
citation frequency In this equation we emphasize that 
the effective parameters can be defined as functionals of 
the self-consistent electron density n(r) as done in the 
Tkatchenko-SchefHer (TS) method p]. Regardless of 
whether one treats a p and uj p as empirical parameters or 
obtains them from the electron density, their values for 
different carbon atoms in polyacene molecules and crys- 
tals turn out to be essentially degenerate (the same holds 
for the hydrogen atoms). This finding can be attributed 
to the rather similar local hybridization environment that 
every atom "feels" inside polyacene molecules. This sim- 
plified model for the polarizability would lead to a similar 
optical absorption spectrum for different polyacenes, in 
stark disagreement with experimental measurements and 
explicit excited-state calculations [H EJ H2] • 

The semilocal approximation for the polarizability in 
Eq. (JlJ neglects the dynamic electric fields that an atom 
experiences from all the other atoms inside a molecule or 
a crystal. Recently, an efficient parameter-free method 
was developed to include these screening effects on the 
polarizability for non-metallic molecules and solids |13) . 
We model the environment as a dipole field and solve 
the resulting classical Dyson-like self-consistent screening 
(SCS) equation [TMH] 
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where a (r;iu>) is the sum of the TS effective atomic 
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polarizabilities [11] . and T(r — r') is the dipole-dipole in- 
teraction tensor (Hartree atomic units are used through- 
out). Equation |2]) is discretized using atomic positions 
as a basis, and then solved directly and exactly by invert- 
ing the tensor corresponding to the coupled dipoles mod- 
eled as quantum harmonic oscillators (QHO). The QHO 
parameters are defined using the TS polarizability |13j . 
The solution of Eq. p| yields the non-local molecular po- 
larizability tensor a^^ (iw), where indices p and q label 
the atoms while indices i and j label the atomic Carte- 
sian coordinates. The contraction of the molecular tensor 
for every atom p yields the atomic polarizability tensors 
a^f(ioj). These tensors now include both the short- 
range hybridization effects from the TS method and the 
long-range response screening from the solution of the 
SCS equation. 

The electrodynamic response included upon solving 
the SCS equation (Eq. [2| allows one to correctly cap- 
ture two important contributions to the polarizability: 
(i) (de)polarization, and (ii) polarizability anisotropy. 
The local "atom-in-a-molecule" polarizability as defined 
in Eq. leads to an essentially isotropic response for 
molecules and solids [13]. The directionality of the 
polarization, well known for polyacenes from experi- 
ments [17] and calculations |12j . emerges from the intrin- 
sic anisotropy of the molecular orbitals and the electro- 
dynamic coupling between them. In the SCS formalism 
of Eq. ([2| , the anisotropy of the molecular polarizability 
stems from the coupling between fluctuating QHOs. For 
a set of small organic molecules, the SCS calculation sig- 
nificantly reduces the error in the molecular anisotropy 
to 23% from 80% in the TS method [13]. The emergence 
of polarizability anisotropy is the main effect brought 
by the inclusion of electrodynamic response effects for 
small gas-phase molecules. In the solid state or for larger 
molecules, the situation is more complex. In a crys- 
talline environment every atom experiences the electric 
field from the atoms within the same molecule (similar 
to the gas phase) , as well as the field produced by neigh- 
boring molecules. While the screening of the molecule 
in the crystal leads to an anisotropic polarizability, one 
also typically finds an appreciable change in the isotropic 
polarizability of the molecule when compared to the gas 
phase. 

In order to assess the relative importance of electrody- 
namic response on the properties of non-covalently bound 
molecular crystals, we have chosen to examine a series 
of polyacene crystals, ranging from naphthalene to pen- 
tacene. This choice allows us to study the evolution of 
response properties of molecular crystals and their sta- 
bility as a function of molecular size and crystal environ- 
ment. Initial crystal structures for each polyacene were 
obtained from the lowest temperature data sets available 
in the Cambridge Structural Database [TS]. The low- 
temperature polymorphs were chosen for tetracene and 
pentacene. The crystal unit cells and the internal geome- 



TABLE I. Unit cell parameters for polyacene molecular crys- 
tals determined from PBE+vdW calculations and low tem- 
perature X-ray experiments. The dielectric constants are re- 
ported using the Clausius-Mosotti equation corresponding to 
Eq. [TUeh), and Eq. [5] (ef u u)- Data are reported for naphtha- 
lene (2 A), anthracene (3A), tetracene (4A), and pentacene 
(5A). 





a (A) b (A) c (A) p (g/ml) 


Efull 


Eh 


2A 


8.fl7 5.897 8.647 


1.244 


3.06 


3.58 


Exp. 


8.f08 5.940 8.647 


1.239 


3.2 [25] 




3A 


8.399 5.906 11.120 


1.313 


3.24 


3.80 


Exp. 


8.414 5.990 11.095 


1.297 


3.2 a 




4A 


6.050 7.706 13.030 


1.343 


3.31 


3.89 


Exp. 6.056 7.838 f3.010 


1.323 






5A 


6.129 7.676 14.531 


1.392 


3.44 


4.08 


Exp. 


6.239 7.636 14.333 


1.397 


2.7-3.89 [mill 


a 



tries were fully optimized using DFT with the generalized 
gradient approximation of Perdew, Burke, and Ernzcrhof 
(PBE) [TS] with vdW interactions treated using the TS 
method [UJ (denoted as PBE+vdW). The CASTEP code 
was used for all calculations [20] [21] . Norm-conserving 
pseudopotentials were employed for carbon, where va- 
lence states included the 2s and 2p electrons. The plane- 
wave basis set cutoff was set to 750 eV, ensuring that 
the total energy and unit cell volume were converged, 
as demonstrated in the study of crystalline indole and 
tetracyanoethylene [22] [23] . The fc-point grid was kept 
to maintain a spacing of 0.07 A -1 . Explicit all-electron 
calculations using the FHI-aims code [53] confirm that 
the binding energies from the pseudopotential calcula- 
tions are converged to better than 0.01 eV per molecule. 

The optimized PBE+vdW lattice parameters for poly- 
acenes along with X-ray measurement results are shown 
in Table HI The overall deviations between our calcula- 
tions and experiments are less than 2% in lattice pa- 
rameters and unit-cell volumes. Similar agreement is 
also found for the internal molecular geometries, where 
the PBE+vdW method predicts the C-C distances with 
an accuracy of 2% in comparison with X-ray measure- 
ments. The polyacene crystal densities predicted by the 
PBE+vdW method are slightly higher than the experi- 
mental ones (except for pentacene), consistent with the 
fact that the inclusion of zero-point energy and thermal 
expansion will decrease the density of the crystal. 

To illustrate the importance of electrodynamic re- 
sponse in polyacene crystals, we have computed the solid 
dielectric constant, e, using the Clausius-Mossotti for- 
mula. The required polarizabilities were obtained from: 
(i) Eq. ([!]), which only includes local hybridization ef- 
fects, and (ii) Eq. which properly accounts for elec- 
trodynamic response screening. Comparison between the 
"hybridized" Eh and the "full" £f u n dielectric constants 
in Table [I] reveals the importance of electrodynamic in- 
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teractions. The £f u n dielectric constants for all polyacene 
solids are decreased by 15% when compared with Eh- The 
fully screened dielectric constants are in excellent agree- 
ment with the measured values of e [U [T71 [251 [SB] . For 
pentacene, £f u n is also close to the values of 3.2 and 3.6, 
obtained by Sharifzadeh et al. [3J from GW calculations 
within the random phase approximation (RPA). Indeed, 
the SCS method solves the RPA equation for a collection 
of QHOs in the dipole approximation. This explains the 
good agreement with GW dielectric constant and also 
the fact that the SCS model successfully reproduces the 
measured dielectric constants of crystalline silicon and 
germanium [I31 [27] . We note that the calculated dielec- 
tric constants from SCS allow us to approximately deter- 
mine the fundamental band-gap of a molecular crystal. 
The fundamental gap of a molecular crystal is reduced 
with respect to that of a molecule in the gas phase, ow- 
ing to the dielectric screening, which reduces the energy 
needed for adding or removing an electron. To evalu- 
ate the fundamental gap of crystalline pentacene we may 
use the polarization model, whereby the gap of the gas 
phase molecule is reduced by 2P. The polarization en- 
ergy, P (in atomic units), is given by — (e — l)/(2i?e), 
where R is determined from the volume per molecule in 
the unit cell as: R — [3V co u/(87r)] 1 / 3 (for two molecules 
per cell) [3J [S81 [25] . Using cell parameters calculated with 
PBE+vdW and £ fu n we obtain P=-1.17 eV. Applying 
the polarization model to the gas-phase gap of 4.57 eV, 
obtained from a GW calculation based on a consistent 
starting point, as described in Ref. [30], we obtain a bulk 
gap of 2.22 eV in good agreement with experiment [3T] 
and with explicit GW calculations for pentacene crys- 
tal [3J. We note that the optical gap is further reduced 
with respect to the fundamental gap due to excitonic ef- 
fects [31 [T21 E] , which are not accounted for by the SCS 
model. 

Experimentally, the stability of molecular crystals is 
measured in terms of their sublimation enthalpy, i.e. the 
energy required to convert a certain amount of molecules 
from the crystalline phase to the gas phase. The sub- 
limation process is carried out at a given temperature 
under constant pressure. The sublimation temperature 
is largely determined by the magnitude of the cohesive 
forces in the crystal. For polyacene crystals, the subli- 
mation temperature varies from w 300 K for naphtalene 
to m 500 K for pentacene [35], illustrating the increase 
in crystal stability for larger acene molecules. There are 
numerous experiments that measure the sublimation en- 
thalpies of polyacenes; in Table |TT] we report a range of 
available values, extrapolated to K [35] • We have only 
taken those values that are recommended as reliable after 
critical revision by the authors of Ref. [32], thus avoid- 
ing anomalously small or large sublimation enthalpies. 
Both naphthalene and anthracene crystals have been vig- 
orously studied, and their sublimation enthalpies are well 
known with a spread of 0.05 eV and 0.12 eV, respectively. 



TABLE II. Lattice energies of polyacene crystals including 
zero-point energy (PBE+MBDh and PBE+MBD calcula- 
tions using optimized PBE+vdW geometries from Table |T]| . 
The range of experimental ("Exp.") "lattice energies" from 
Ref. [32) extrapolated to K. Also shown are the AH from 
room temperature (RT) to K calculated from an integral 
over experimental C P (T) data, and the harmonic zero-point 
energy calculated using the PBE+vdW method. All values 



are 


in units of eV. 












A#(RT-> 0) 


ZPE 


PBE+MBDh PBE+MBD 


Exp. 




2A 


0.041 


0.069 


-0.993 


-0.862 


-0.803 to - 


-0.752 


3A 


0.034 


0.078 


-1.433 


-1.206 


-1.148 to - 


-1.024 


4A 


-0.016 


0.110 


-1.951 


-1.587 


-1.525 to - 


-1.299 


5A 


-0.032 


0.115 


-2.501 


-2.018 


-2.082 to - 


-1.533 



There are fewer measurements available for tetracene and 
pentacene, and for the latter the three available experi- 
mental values deviate by 0.55 eV. 

In order to compare theoretical lattice energies to the 
measured sublimation enthalpies, the sublimation en- 
thalpies need to be extrapolated to zero temperature 
by adding the enthalpy difference [i?°(T) - ff°(0)] - 
[Hg(T) — Hg(0)], where the subscript c refers to the 
crystal, whereas g refers to the gas phase. In this work, 
we have calculated this enthalpy difference by integrat- 
ing the measured heat capacity C P (T) for acene crys- 
tals, and the extrapolated gas-phase heat capacity from 
Refs. [33H36] . Even at K, the sublimation enthalpy in- 
cludes zero-point vibrational effects, and these have to 
be considered when comparing calculated lattice ener- 
gies to the experimental enthalpy extrapolated to K. 
Here we determined the zero-point energy (ZPE) from 
phonon calculations using the PBE+vdW method with 
the supercell formalism in CASTEP [20]. Special care 
has been taken to converge the supercell size and plane- 
wave cutoff when performing phonon calculations. We 
estimate that our ZPE calculations are converged to 5 
meV/molecule. Note that vdW interactions contribute 
significantly to the ZPE energy and have to be included 
in phonon calculations to reach this level of accuracy. 

Now we analyze the impact of electrodynamic re- 
sponse on the sublimation enthalpies of polyacene crys- 
tals. For this purpose we combine the PBE func- 
tional with the recently developed many-body dispersion 
(MBD) method [T3J. The MBD energy expression com- 
putes the long-range many-body dispersion energy to in- 
finite order for molecules and solids with a finite band 
gap. In the MBD method, the full electronic system is 
mapped to a system of quantum harmonic oscillators. 
The QHO polarizabilities can be obtained either from 
Eq. (from now on called the DFT+MBDh method) 
or Eq. @ (from now on called the DFT+MBD method), 
thus allowing us to clearly assess the effect of electrody- 
namic screening on the dispersion energy and stability 
of molecular crystals. The ZPE- inclusive lattice energies 
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FIG. 1. (Color online) Left: The difference in the vdW energy 
between the PBE+MBD method including electrodynamic 
response and the PBE+MBDh method based on "atom-in- 
a-molecule" polarizabilities of Eq. |TJ. Right: The difference 
in the C 6 coefficients between PBE+MBD and PBE+MBDh 
methods per carbon atom. 



obtained with both methods are shown in Tabic |TT] and 
compared with experimental data. 

The inspection of Table [TT| illustrates the crucial impor- 
tance of electrodynamic response for the stability of poly- 
acene crystals. The PBE+MBDh method uses semilocal 
hybridized polarizabilities and overestimates the experi- 
mental sublimation enthalpies by more than 0.20 eV for 
naphthalene and up to 0.44 eV for pentacene. Upon 
including the response screening, as depolarization re- 
duces the stability of naphthalene by 0.13 eV and of 
pentacene by 0.48 eV compared to the PBE+MBDh ap- 
proach. To assess the influence of the underlying DFT 
functional, we have also carried out calculations using 
the PBE-based hybrid functional, PBE0 [33135], com- 
bined with the MBD method. The PBE0 functional de- 
scribes the electrostatic and inductive intermolecular in- 
teractions more accurately |13[ 139]. However, we found 
that the PBE0+MBD approach yields essentially the 
same lattice energies as the PBE+MBD method for all 
polyacenes (within 10 meV per molecule). 

The remaining slight overestimation of lattice energies 
in Table [IT] compared to the experimental range can be 
explained by the fact that the sublimation enthalpy is 
measured at finite temperature, where the unit cell un- 
dergoes thermal expansion. When using the experimen- 
tal unit cell at 295 K for naphthalene, the PBE+MBD 
method yields a lattice energy that is increased by 50 
meV, essentially within the experimental range reported 
in Table HU 

Finally, we explain the observed difference between 
the PBE+MBDh and PBE+MBD methods by analyz- 
ing the vdW dispersion energies in the gas and crystal 
phases in Figure [T] Along with vdW energies, we also 



show the change in the vdW C@ coefficient per carbon 
atom upon inclusion of the electrodynamic response in 
the MBD method. For the gas phase molecules, there 
is a significant dipole polarization along the long molec- 
ular axis, which increases the molecular polarizabilities 
and Cq coefficients, leading to an increase in the vdW 
energy in the PBE+MBD method when compared to 
PBE+MBDh. However, since the vdW energy contri- 
bution to the molecular stability is relatively small, the 
change due to electrodynamic response is only -0.05 eV 
for naphthalene and up to -0.13 eV for pentacene. The 
electrodynamic response gives rise to a radically differ- 
ent situation in the crystal phase; in this case the in- 
teraction with neighboring molecules leads to overall de- 
polarization, decreasing the Cg coefficients by roughly 
a constant amount, when compared to the gas-phase 
molecules. However, the vdW energy makes a larger 
absolute contribution to the stability of acene crystals. 
This explains the sharp decrease of the vdW energy in 
the crystal predicted by the PBE+MBD method when 
compared to the PBE+MBDh approach. Overall, the 
opposite effect of screening for the gas-phase molecule 
and the crystal explains the significant reduction of the 
lattice energy, shown in Table [Hj upon including electro- 
dynamic response in the PBE+MBD method. 

In summary, we have quantitatively established a con- 
nection between collective electrodynamic response and 
the stability of molecular crystals. Our results demon- 
strate that molecular crystals are significantly more com- 
plex than a simple collection of constituent molecules, 
and that electrodynamic response is crucial for explain- 
ing many of the unique properties of molecular crystals. 
We provide to our knowledge the first quantification of 
the influence of electrodynamic response on sublimation 
enthalpies of polyacene crystals, and it is not unreason- 
able to expect that our findings will hold in general for 
other classes of molecular solids. 

A.T. acknowledges support from the European Re- 
search Council (ERC Starting Grant VDW-CMAT). B.S. ac- 
knowledges support from the Eberly Science Foundation. 



[1] 

[2] 



[3] 



tkatchenko@fhi-berlin.mpg.de 
S. R. Forrest, Nature 428, 911 (2004) 
T. W. Kelley, P. F. Baude, C. Gerlach, D. E. Ender, 
D. Muyres, M. A. Haase, D. E. Vogel, and S. D. Theiss, 
Chem. Mater. 16, 4413 (2004) 

S. Sharifzadeh, A. Biller, L. Kronik, and J. B. Neaton, 
Phys. Rev. B 85, 125307 (2012) 
[4] M. Pope and C. E. Swenberg, Electronic Processes in Or- 
ganic Crystals and Polymers (Oxford University Press, 
New York, 1999) 
[5] J. E. Anthony, Angew. Chem. Int. Ed. 47, 452 (2008) 
[6] J. S. Toll, Phys. Rev. B 104, 1760 (1956) 
[7] M. A. Neumann and M. A. Perrin, J. Phys. Chem. B 
109, 15531 (2005) 



5 



[8] D. A. Bardwell et ai, Acta. Cryst. B 67, 535 (2011) 
[9] S. Wen, K. Nanda, Y. Huang, and G. Beran, Phys. Chem. 
Chem. Phys. 14, 7578 (2012) 
[10] A. J. Stone, The theory of intermolecular forces (Oxford 

University Press, 1996) 
[11] A. Tkatchenko and M. Scheffler, Phys. Rev. Lett. 102, 
073005 (2009) 

[12] M. L. Tiago, J. E. Northrup, and S. G. Louie, Phys. Rev. 
B 67, 115212 (2003) 

[13] A. Tkatchenko, R. A. DiStasio, Jr., R. Car, and M. Schef- 
fler, Phys. Rev. Lett. 108, 236402 (2012) 

[14] B. U. Felderhof, Physica 76, 486 (1974) 

[15] D. W. Oxtoby and W. M. Gelbart, Mol. Phys. 29, 1569 
(1975) 

[16] B. T. Thole, Chem. Phys. 59, 341 (1981) 

[17] D. Faltermeier, B. Gompf, M. Dressel, A. K. Tripathi, 

and J. Pflaum, Phys. Rev. B 74, 125416 (2006) 
[18] F. H. Allen, Acta. Crys. B58, 380 (2002) 
[19] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 

Lett. 77, 3865 (1996) 
[20] S. J. Clark, M. D. Segall, C. J. Pickard, P. J. Hasnip, 

M. J. Probert, K. Refson, and M. C. Payne, Zeit. Krist. 

220, 567 (2005) 
[21] E. R. McNellis, J. Meyer, and K. Reuter, Phys. Rev. B 

80, 205414 (2009) 
[22] B. Schatschneider and J. -J. Liang, J. Chem. Phys. 135, 

164508 (2011) 

[23] B. Schatschneider, J.-J. Liang, S. Jezowksi, and 
A. Tkatchenko, CrystEngComm 14, 4656 (2012) 

[24] V. Blum, R. Gehrke, F. Hanke, P. Havu, V. Havu, 
X. Ren, K. Reuter, and M. Scheffler, Comp. Phys. Comm. 
180, 2175 (2009) 



[25] T. Suthana, N. Rajesha, P. Dhanaraja, and C. Mahade- 
van, Spectrochimica Acta Part A 75, 69 (2010) 

[26] C. H. Kim, O. Yaghmazadeh, D. Tondelier, Y. B. Jeong, 
Y. Bonnassieux, and G. Horowitz, J. Appl. Phys. 109, 
083710 (2011) 

[27] V. V. Gobre et al, to be published 

[28] I. G. Hill, A. Kahn, Z. G. Soos, and R. A. Pascal, Jr. 

327, 181 (2000) 
[29] J. B. Neaton, M. S. Hybertsen, and S. G. Louie, Phys. 

Rev. Lett. 97, 216405 (2006) 
[30] T. Korzdorfer and N. Marom, Phys. Rev. B 86, 041110 

(2012) 

[31] E. A. Silinsh, V. A. Kolcsnikov, I. J. Muzikante, and 

D. R. Balode, physica status solidi (b) 113, 379 (1982) 
[32] M. V. Roux, M. Temprado, J. S. Chickos, and Y. Nagano, 

J. Phys. Chem. Ref. Data 37, 1855 (2008) 
[33] NIST WebBook and references within, 

http://webbook.nist.gov/chemistry/name-ser.html. 
[34] J. P. McCullough, H. L. Finke, J. F. Messerly, T. C. 

Kincheloe, and G. Waddington, J. Phys. Chem. 61, 1105 

(1957) 

[35] P. Goursot, H. L. Girdhar, and E. F. Wcstrum, Jr., J. 

Phys. Chem. 74, 2538 (1970) 
[36] W.-K. Wong and E. F. Westrum, Jr., Mol. Crys. Liq. 

Cryst. 61, 207 (1980) 
[37] J. P. Perdew, M. Ernzerhof, and K. Burke, J. Comp. 

Phys. 105, 9982 (1996) 
[38] C. Adamo and V. Barone, J. Comp. Phys. 110, 6158 

(1999) 

[39] B. Santra, J. Klimes, D. Alfe, A. Tkatchenko, B. Slater, 
A. Michaelides, R. Car, and M. Scheffler, Phys. Rev. Lett. 
107, 185701 (2011) 



